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ABSTRACT 

Aims. Variability at all observed wavelengths is a distinctive property of active galactic nuclei (AGN). Hard X-rays provide us with a 
view of the innermost regions of AGN, mostly unbiased by absorption along the line of sight. Characterizing the intrinsic hard X-ray 
variability of a large AGN sample and comparing it to the results obtained at lower X-ray energies can significantly contribute to our 
understanding of the mechanisms underlying the high-energy radiation. 

Methods. SwiftfBAA provides us with the unique opportunity to follow, on time scales of days to years and with a regular sampling, 
the 14—195 keV emission of the largest AGN sample available up to date for this kind of investigation. As a continuation of an early 
work on the first 9 months of BAT data, we study the amplitude of the variations, and their dependence on sub-class and on energy, 
for a sample of 1 10 radio quiet and radio loud AGN selected from the BAT 58-month survey. 

Results. About 80% of the AGN in the sample are found to exhibit significant variability on months to years time scales, radio loud 
sources being the most variable, and Seyfert 1.5-2 galaxies being slightly more variable than Seyfert 1, while absorbed and unabsorbed 
objects show similar timing properties. The amplitude of the variations and their energy dependence are incompatible with variability 
being driven at hard X-rays by changes of the absorption column density. In general, the variations in the 14-24 and 35-100 keV bands 
are well correlated, suggesting a common origin of the variability across the BAT energy band. However, radio quiet AGN display on 
average 10% larger variations at 14-24 keV than at 35-100 keV, and a softer-when-brighter behavior for most of the Seyfert galaxies 
with detectable spectral variability on month time scale. In addition, sources with harder spectra are found to be more variable than 
softer ones, opposite to what it is observed below 10 keV. These properties are generally consistent with a variable, in flux and shape, 
power law continuum, pivoting at energies > 50keV, to which a constant reflection component is superposed. When the same time 
scales are considered, the timing properties of AGN at hard X-rays are comparable to those at lower energies, with at least some of 
the differences possibly ascribable to components contributing differently in the two energy domains (e.g., reflection, absorption). 
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1. Introduction 

Besides mapping the innermost regions of active galactic nuclei 
(AGN), the hardest X-rays >15 keV present the further advan- 
tage of not being affected by absorption, provided that the hydro- 
gen column density Ah along the line of sight is lower than a few 
10 23 cm -2 . Hence, hard X-rays can be effectively used to study 
the intrinsic properties of AGN emission, testing the validity of 
unification models that in first approximation explain the differ- 
ences between the different AGN classes as a function of the 
viewing angle (e.g.. lBeckmann & Shradeili201 2i) . In radio quiet 
AGN, the hard X-ray emission is postulated to originate in a two 
phase medium, where the soft disk photons are C omp tonized 
by ele ctrons in a hot plasma above the disk ( Ha ardt & Maraschil 
119931) . In radio loud objects, an additional contribution from in- 
verse Compton scattering from relativistic electrons in the jet as 
well as synchrotron radiation can dominate the X-ray output. 


Moreover, AGN are known to exhibit variability at all ob- 
served frequencies, whose study can provide important infor- 
mation about the physics, the structure and the dynamics of 
the emitting source. In the softer X-ray domain (< lOkeV), 
extensive variability studies have been carried out, especially 
in the recent years thanks to the numerous monitoring cam- 
paigns with RXTE and XMM-Newton , covering from hour 
to year time scales. Many important results have been ob- 
tained with these observations, as for instance the complex 
correlation o f the X-rays with the emission at other wave- 
lengths (e.g., McHardv et al. 2007t Soldi et al. 20081: iBell et alJ 


1201 lb iMehdipour et al. 201 lb Chatteriee et al. 20JJ]T in par- 
ticular in the optical-UV band (e.g., Shemmer et al. 2003; 


IUttlevllT005b lArevalo et al.ll2008bb Breedt et al.ll2010l): the study 


of X-ray lag spectra (e.g., lArevalo et al. 2008ab IZoghbi et all 
12010b iDe Marco et alJl201 3h : the correlatio ns between variabil- 
ity properties and black ho l e mas s (e.g., Czerny et alJ 1200 lb 


iLu & Yul 1200 lb llJttlev et al 
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cretion rate ( McHardv et alJ 


Panadakis 20041) ar id ac- 
Kording et al. 200 7). See 
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iMcHardvl (1201 Ol) for a review on X-ray variability and its scal- 
ing from stellar to super-massive black holes. On the other 
hand, the large majority of the hard X-ray variability studies 
performed in the last years focused on the spectral variability 
of AGN at different flux levels for single bright sources or for 
small samples of selected o bjects, thanks to pointed observa- 


2000; 

de Rosa et al. 2007J), INTEGRAL (e.g., Pian et al. 2006, 

2011; 

Beckmann et al. 2008; Lubiriski et al. 2010‘ Soldi et al. 

2011a 

Petrucci et al. 20131) and Suzaku (e.g.. Reeves et al. 2007; 

Itoh et al.l! 

2008 

Terashima et al.l 2009; iFukazawa et al. 2011; 

Reis et al. 

2012 

) satellites. These X-ray studies suggest differ- 


ent origins for the observed variability, as for example changes 
of the amount of absorption or of its ionisation, fluctuations in 
the seed photon flux, intrinsic modifications of the corona prop- 
erties or geometry (or of the jet parameters, in case of radio loud 
AGN), a variable contribution of different spectral components, 
with at least some of these scenarios being driven by accretion 
rate variations. 


However, with the only exception being the CGRO/BATSE 
instrument, detecting only a handful of AGN (lHarmon et al.l 
2004 ), hard X-ray, long-term monitoring has not been possi- 
ble until recent years, because of the observing strategy of the 
hard X-ray satellites and the limited size of the field of view 
of their instruments. Furthermore, last generation and future 
hard-X-ray telescopes such as NuSTAR will have small fields 
of view and relatively inflexible scheduling capabilities, so fu- 
ture studies may tend to be limited to either short time scales 
or in the number of observed objects. Since November 2004, 
the Burst Alert Telescope (BAT, |Barthehriv_e^if[ 20051) instru- 
ment on board the Swift satellite dGehrels et alJ 20041) has been 
observing the sky in the 14-195 keV energy range. Thanks to 
its large field of view of ~ 1.4 sr and to Swiff s observing strat- 
egy, the BAT has been monitoring a large number of hard X-ray 
sources (Cusuman o et aklfeOlOt Baumgart ner et al. 11201 3l) . pro- 
viding for the first time a long-term and sufficiently sampled 
data set. Swift/BAT light curves on different time scales, energy 
bands and extracted with different methods are provided through 
the hard X-ray transient monitoring pages[] (iKrimin et al.ll2013l) . 
and for the BAT 5 8 -month hard X-ray survejQ including more 
than 1000 objects, about 60% of which are AGN. The latter 
has been recently update d to cover 70 months of observations 
dBaumgartner et al.ll2013l) . 


A first study of the variability of hard X-ray selected AGN 
using BAT data has been presented bv lBeckmann et all ( 2007 a). 
It included 44 AGN detected with high significance over the first 
9 months of BAT observations. The variability study in the 14- 
195 keV band has been performed using a maximum likelihood 
estimator and the structure function analysis. Among the main 
results of this work, Seyfert 2/obscured objects were found to 
be more variable than Seyfert 1/unobscured ones, and an anti- 
correlation has been detected between variability and luminos- 
ity. 

Preliminary studies of the 5 years BAT observations but with 
limited numbers o f AGN have been presented by ISoldi et al.l 
d2010l) . iRicci et akl (1201 lal) and ICaballero-Garcia et al.l (120121) . 
In particular the latter study focused on 5 bright AGN, detecting 
spectral variability in three of them, consistent with the electron 
plasma temperature decreasing with increasing flux, within the 
Comptonization scenario. 


1 http://swift.gsfc.nasa.gov/docs/swift/results/transients 

2 http .'//swift, gsfc. nasa. gov/docs/swift/results/bs58mon/ 


IShimizu & Mushotzkv ( 2013L hereafter SM13) have calcu- 
lated for the first time the power density spectra (PDS) of 30 
AGN at hard X-rays (all these objects are also included in our 
study), using BAT 14-150 keV data covering time scales from 
one week to a few years (vpds ~ 1 0 x - 10 6 Hz). All but 
one PDS were found to be well fitted by a simple power law 
with slope a ~ -0.8, similarly to the PDS measured in the 
2-10 keV domain at time frequencies lower than a break fre- 
quency 1/T B . In fact, the PDS of AGN in the 2-10 keV do- 
main has a characteristic power law shape with index apos = 
-1 for long time scales, steepening to apos = -2 for times 
scales shorter than a source- spe cific break time 'If ( UttlevlfT005i. 
iGonzalez-Martm & Vaughan»2012h . which is found to scale with 
black hole mass and accretion rate following the relation de- 
scribed by iMcHardv et al.l d2006l) . The lack of correlations be- 
tween variability and luminosity, and variability and black hole 
mass in the sample presented by IShimizu & Mushotzkvi (1201 3l) 
has, therefore, been ascribed to the long time scales probed by 
this study. No decisive evidence was found to determine whether 
or not AGN long-term variability is energy dependent. 

We present here our variability analysis of the largest AGN 
sample with a sufficient hard X-ray monitoring to perform this 
kind of study. Our anaiysis provides an orthogonai approach 
to that presented by SMI 3, as we consider the frequency- 
independent variations of the full long-term light curves, and 
therefore we are able to extend our study to a much larger AGN 
sample, while for the PDS analysis a high signal-to-noise ratio 
is required, thus limiting the use of this technique to the 30 ob- 
jects already presented by SM13. Furthermore, we are also able 
to investigate the dependence of variability on energy in the hard 
X-ray band. 

A description of the AGN light curves in the BAT 58-month 
survey is presented in Sect. [2] together with the estimate of the 
systematic uncertainties. In Sect.[3]the maximum likelihood esti- 
mator for variability is introduced and the selected sample of 1 10 
AGN is described. The correlation analysis between variabil- 
ity and the AGN fundamental parameters is detailed in Sect. |4] 
while Sect. [5]presents the study of spectral variability within the 
BAT energy range. We discuss our results (Sect. [6]) with partic- 
ular emphasis on the comparison between hard X-ray and soft 
X-ray variability, since this can provide important insights on 
the mechanisms responsible for the observed emission across the 
high-energy spectrum. A summary of our conclusions is given in 
Sect. |7] 

Considering the importance of the analysed time scales in the 
results presented here, we will generally refer to long time scales 
(and low time frequencies) for variability studies considered to 
sample the PDS above the break 7 b, and to short time scales 
(and high time frequencies) below the break. In the following 
we use a ACDM cosmology with Om = 0.3, Qa = 0.7 and 
Ho = 73kms~ 1 2 Mpc~ 1 . 


2. Swiff/ BAT light curves analysis 

Among the 1092 sources detected in the BAT 58-month sur- 
vey, 625 are AGN. We exclude from this list the 6 pairs of 
AGN whose sky positions are too close to safely avoid con- 
tamination of the light curves, i.e., objects with BAT coordi- 
nates within 2 arcmin from each other. For these 613 AGN we 
retrieved the publicly available light curves covering the time 
from the beginning of the mission up to May 2010, therefore in- 


Article number, page 2 ofQT] 


S. Soldi et al.: Long-term variability of AGN at hard X-rays 


eluding up to 66 months of dat;Q. The light curves have been 
extracted from the snapshot (i.e. a single Swift pointing lasting 
about 20 minutes) images corrected for off-axis effects. The im- 
ages were created in 8 energy bands, i.e. 14-20, 20-24, 24-35, 
35-50, 50-75, 75-100, 100-150, 150-195 keV. In order to con- 
vert count rates in the total 14-195 keV band to fluxes, we as- 
sume a Crab-like spectrum (T = 2.15: lTueller et al .1120 Hit) with a 
flux of fo-ab, 14 — 195 keV = 2.44 x 10 -8 erg cm -2 s - 1 corresponding 
to a count rate of xc m b,i 4 -i 95 kev = 0.0418 counts s -1 detector -1 
(average from the Crab snapshot light curve). The complete de- 
scription o f the BAT dat a analy sis and light curve extraction is 
provided in lTueller et al.l (1201 Pi) and lBaumeartner et al.l (1201 3h . 

The reference light curves used for our analysis are obtained 
from the combined 14-195 keV light curves first rebinned to 
30 days and then filtered to exclude points with exposure time 
shorter than one day, and data points with error bars larger than 
AF cutj 30 days =1.7 mCrab. To derive this value, the histogram of 
the flux uncertainties for all light curves was built in logarith- 
mic space, and the histogram peak and minimum identified. The 
value AT cut , 30 days = 2 x ATp e ak - AC m j n (i.e., as far above the 
peak of the distribution as the minimum is below) is considered 
to mark the beginning of the high-value tail of the distribution 
of flux uncertainties and, therefore, is used as the cut-off value, 

AT C ut,30days ■ 

2.1. Light curve systematic errors 

Residual systematic uncertainties can affect the light curves and 
need to be taken into account before a variability analysis can 
be performed. Among them, the most important contributions 
come from a component cr sys a = sysA x F which is proportional 
to the source flux, a component cr sys B = .vy,v H x cr sta t that can be 
derived from empty sky positions, and a term cr sys c = sysc x 
cr stat due to summing up energy and time bins that are somehow 
correlated. cr sys a is estimated using the Crab light curve under 
the assumption that it is constant. Fitting the Crab orbital light 
curve with a constant, a cr sys A = 0.06 x F systematic error needs 
to be added in order to obtain a reduced x~ equal to 1 . A further 
source of uncertainty is expected following the finding of long- 
term variability in the Crab lig ht curve (up to ~7% on year time 
scale: lwilson-Hodge et alJl201 ill , therefore, the value of the sys& 
term is likely to be overestimated here. Due to the difficulty of 
disentangling the systematic component from the intrinsic Crab 
variability, we prefer to choose a conservative approach and to 
adopt the sysA = 0.06 value. In any case, even when reducing 
this term to 0.01, the amplitude of the variations increases on 
average by only 2-3%, largely within the estimated uncertainties 
on the variability estimator (see Sect. ED . 

The component cr sys B is estimated by constructing orbital 
light curves from blank sky positions and building the histogram 
of all the resulting significances. By selecting sky positions with 
no sources, a Gaussian distribution centered on zero for the sig- 
nificances is expected. Any significant deviation from a zero 
mean value or from a width equal to 1 would point towards sys- 
tematic uncertainties. Indeed the statistical errors need to be in- 
creased by cr sys B = 0.08 x cr stat in order to obtain a Gaussian fit 
to the histogram of significances with width equal to 1. 

Finally, from the comparison between the 14-195 keV 
monthly light curves extracted from the 1 -month integrated mo- 
saic images and the monthly light curves obtained by rebin- 

3 We note that even though 58 months of BAT data were used for 
source detection, the light curves used in this work extend beyond this 
data set, i.e. up to 66 months. 
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Fig. 1. Examples of 30-d binned, 14-195 keV BAT light curves for ob- 
jects presenting low to high variability (from top to bottom) and with 
different signal-to-noise ratios. The estimator S v represents the intrin- 
sic variability of a light curve, renormalized to the average flux and in 
percentage units (see Sect. 13. il l. 


ning the orbital light curves and summing up 8 energy bands 
over the full 14-195 keV range, we estimate that an additional 
statistical error contribution, cr sys c = 0.08 x cr stat , is needed 
to take into account the effects of correlated energy and time 
bins. The final statistical plus systematic error is given by <x tot = 
^J(sys A x F) 2 + [cr stat ( 1 + sy.s B + vysc)] 2 - 

In Fig.Q]five example light curves are presented with differ- 
ent signal-to-noise ratios and variability. 


3. Hard X-ray variability properties 

As an initial estimate of the variability of our AGN sample, we 
fit the 30-day binned light curves of the 613 AGN with a constant 
function, and apply a^ 2 test. The variability index V is defined as 
V = - log( 1 -P x - ’■), where P x i is the null hypothesis probability to 
obtain such a^ 2 if the source were constant. For 36% of the light 
curves in the sample the fit t o a constant intensity source results 
in V > 1.3 (i.e., P x i < 5%; Pao lillo et alj|2004t lLanzuisi et al.l 
1201 3l) . indicating that these objects exhibit significant variability 
(Fig. [2]). In particular, 34% of the Seyfert galaxies are found to 
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Fig. 2. Histograms of the time-averaged 14-195 keV flux (top), of the 
variability index V from the^ 2 test (middle; the vertical dashed line in- 
dicates V = 1.3), and of the variability amplitude estimator S v (bottom) 
for the total sample of 613 AGN (black) and for the selected sample of 
110 AGN (red hatched). We assigned V = 10 to those objects with 
P x 2 < 1(T 6 . S v could be computed for 464 of the BAT detected AGN. 


be variable against the 44% of the radio loud population. Among 
Seyfert galaxies, about 41% of type-2 objects show variability, 
contrary to only 29% of the type-1 objects. 


3.1. Amplitude of the variations 

Several methods can be used to quantify the amplitude of the 
variations in a light c urve, e.g., excess var i ance or fractional 
variability ampli tude (IVaughan et al.l 12003b iPonti et all l2004t 
ISoldi et akll2008l) . lAlmaini et al. ~ (l2000l) proposed a numerical 
approach arguing that this is a more appropriate method in pres- 
ence of non uniform measurement uncertainties, while recently 
lAllevato et al.l (1201 3l) showed that it provides results equivalent 
to the normalized excess variance. The Almaini method is based 
on a maximum-likelihood estimate of the ctq parameter, repre- 
senting the variability of the light curve and which, in the case of 
constant measurement uncertainties (<x, = constant), reduces to 
the excess variance ( <xq = cr xs ). The crq_parameter is defined so 
as to satisfy the followi ng equation (see Almaini et al. 200(3 and 
iBeckmann et al.ll2007al for more details): 


y, [U'i - X) 2 ~ (Of + 0 "q)] 

h (-r + ^) 2 


(i) 


where x is the mean value of the light curve x, with mea- 
surement uncertainties it, . The maximum-likelihood estimator 
has b een successfu l ly use d also in, e.g., [Mateos et al. (2007) 
and in Chitnis et af] ( 20091) . For continuity and comparison with 
IBeckmann et al.l ( 2007al) . we choose to apply here the Almaini 
method. Unlike the original prescription, we use the weighted 


rather than the arithmetic mean to compute ctq, because it pro- 
vides in general more stable results, being less sensitive to the 
presence of outlier points. Yet, we verified that equivalent re- 
sults are obtained in the two cases when analysing the final AGN 
sample presented here. 

In order to compare the variability of different objects, we 
renormalize cr q to the average flux (F) of the source, obtaining 
the Sy = <xq /(F) x 100% variability estimator (Fig. [2}- There- 
fore, S v measures the amplitude of the intrinsic variability of a 
light curve, corrected by measurement uncertainties, renormal- 
ized to the average flux and in percentage units. Differently from 
IBeckmann et alJ (12007 ah . the systematic uncertainties deduced 
from the Crab and blank sky positions have already been in- 
cluded in the orbital light curves and therefore, no further cor- 
rection to S v is required. 

In computing S y the uncertaint y (X, n ,-i S is determined with 
a bootstrap technique (ISimpson & Maver-Hasselwanderil 1 986i) . 
For each object, 100,000 light curves are randomly drawn from 
the original one, with the same total number of points and al- 
lowing the same data point to be drawn more than once. The 
variability estimator S v is then computed for each simulated 
light curve and its frequency distribution built. The 15.9th and 
84.1th percentiles of this distribution are taken as an estimate of 
the lcr confidence intervals (see also lVaughan et al1l2003l) . Due 
to the stochastic nature of variability, a further source of uncer- 
tainty cr samp i is introduced by the uneven sampling of the light 
curves. This contribution can be understood as the dispersion 
of .S' v values measured on the same light curve on which dif- 
ferent samplings are applied. We estimate this contribution with 
the bootstrap method but this time using the initial, unfiltered 
light curves (i.e., all having regular 30-day sampling). For each 
object, 10,000 light curves with N points are randomly drawn 
from the original one, where N is the number of points of the 
filtered light curve, and where every data point cannot be drawn 
more than once. The 15.9th and 84. 1th percentiles of the S y fre- 
quency distribution are taken as an estimate of the lcr confidence 
intervals. The final uncertainty on S y is obtained by combining 
in quadrature cr meas and cr samp i. 

3.2. Sample selection 

The final sample of 110 AGN is selected based on the follow- 
ing two criteria: 1) the average value of the signal-to-noise ratio 
of the points in the light curve is larger than 2; 2) Sy can be 
measured for the given 30-day binned, filtered light curve, i.e., 
the uncertainties are smaller than the measurable intrinsic vari- 
ability. The first criterion is essentially (but not exactly) equiv- 
alent to using a flux threshold of F( 14 - 195 keV) > 3.8 x 
10 _n erg cm -2 s -1 or a signal-to-noise ratio over the full light 
curve larger than 15, and it selects 1 15 sources (Fig. 0 . The sec- 
ond criterion excludes only 5 additional sources. The filtering 
applied to the light curves (see Sect. 0 excludes no more than 
22% (~10% on average) of the data for the objects in the se- 
lected sampleQ The resulting sample is composed of 88 Seyfert 
galaxies (36 Seyfert 1, 17 of intermediate type, 32 Seyfert 2, and 

3 Narrow Line Seyfert 1), 9 radio galaxies and 13 blazars. These 
sources cover a range of 14-195 keV fluxes between F = 10 -11 

4 The values of 5 v computed on the filtered and unfiltered light curves 
are well compatible with an average scatter of < 2% for the selected 
sample, well within the estimated uncertainties on Sv- On the other 
hand, for the BAT AGN excluded from the analysis the scatter is on 
average 4 times larger, further justifying our choice to limit the analysis 
to a sample of light curves for which the variability measurement is not 
so strongly influenced by single data points. 
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and 10 erg cm 2 s 1 (Fig. O and redshifts up to z = 2.5 (see 
Fig. 1 in lSoldi et alj201 j for the redshift distribution of the sam- 
ple). The non-blazar objects belong to the local AGN population 
with redshifts up to z = 0.1 and an average value of (z) = 0.024 
(Table Q}. 

In Table [2] the values of S v are reported for the 110 AGN 
selected for the variability amplitude analysis. The blazar pop- 
ulation shows a larger average variability ((Sv)bia = 33 + 2%) 
compared to the radio quiet objects ((Sv)sey = 19.3 + 0.5%). A 
Kolmogorov-Smirnov test provides a probability Pks 4 98% 
that the blazar and Seyfert samples are not drawn from the 
same parent population. Yet, extreme cases of high variability 
are present in both the blazar and Seyfert classes. An average 
variability S v > 90% is in fact detected in the two gamma-ray 
bright blazars QSO B0033+595 andMrk421 and in the Seyfert 2 
galaxy 2MASX J04440903 +28 13003. 

The radio galaxies present an intermediate behavior between 
blazars and Seyferts, with (.S' v)rg = 24.0 ± 1.4%. Within the 
Seyfert class, intermediate and type 2 objects ((Sv)seyi. 5+2 = 
20.8+0.6%) are found to be only slightly more variable than type 

1 ((Sv)seyi = 16.5+0.9%), and a KS-test probability P K s ^ 94% 

indicates that the Seyfert 1 and Seyfert 1 .5-2 samples could be 
drawn from the same parent distribution. Furthermore, no dif- 
ference is found in the average variability of absorbed and un- 
absorbed Seyfert galaxies (Ah = 10 22 cm~ 2 being the dividing 
line), with (S v ) abs = 19.7 + 0.6% and (Sv)unabs = 19 ± 1%. 

The average properties of the different AGN classes are sum- 
marized in Table Q] 


3.3. Black hole mass , hydrogen column density and 
Eddington ratio 

In order to test the dependence of variability on different AGN 
parameters, we collected from the literature the values of black 
hole masses Mbh and hydrogen column density Ah for our sam- 
ple. The masses have been calculated with different methods, 
and when not available in the literature, the uncertainties on 
the values were estimated from the typical accuracy provided 
by the method following iBeckmann et alJ ( 2009 ). Out of the 
110 objects, 92 have mass estimates (only an upper limit for 
MR 2251-178; Table 0, 86 have Ah values, and for 23 only 
Ah upper limits could be found. 

For 63 objects (56 Seyferts, 4 radio galaxies and 3 blazars) 
we collected the bolometric luminosities from the literature, se- 
lecting tho se estimated from the fitting of the spectral energy dis- 
tributions (IWoo & Urrv||2002l: IVas udevan & Fabianll2007l 12009: 


IVasudevan et alJl2009[ 201 dt Table [21. In order to estimate the 
bolom etric luminosity from the X-ray one, co nstant correc- 
tions dWinter et all 120091: IBeckmann etak 2009), and correc- 
tions depending on the luminosity (lM arconi e t all 120041) . on 
the Edding ton ratio ( Vasudevan & Fab ian 2007), on the spec- 
tral shape dFabian et al.1 2009 ). or on the intrinsic absorption 
d Vasudevan et 'al J 120 1 Ol) have been proposed, but all fail to re- 
produce the large scatter of bolometric versus X-ray luminosity. 
However, in order to extend the sample for which we are able to 
compute the Eddington ratio, we also considered two of those 
bolometric corrections. The first one follows IBeckmann et al .1 
(120091) - i.e., assumes that the bolometric luminosity L ho | is 2 
times larger than the 1-200 keV luminosity, where L\ 200 kev = 
1.8 x Li 4 _i 95 kev for a simple power law spectrum with photon 
index F = 1.9. The second one is based on the relation provided 
by iMarconi et al.1 d2004l) to compute Lboi using a luminosity- 
dependent bolometric correction. Since these corrections allow 
to have estimates of L ho | for additional 29 objects (19 Seyfert) 



Fig. 3. Variability estimator Sy versus luminosity in the 14-195 keV 
band for the Seyfert galaxies in our sample. Yellow triangles represent 
the average of the variability estimator in different luminosity bins fol- 
lowing j^eyatoeFaT,: (12013). 



Fig. 4. Variability estimator Sy versus intrinsic absorption for the 
Seyfert galaxy sample. Upper limits on Ah are indicated with arrows. 
Yellow triangles represent the average of the variability estimator in dif- 
ferent absorption bins. 


and they trace sufficiently well the bolometric luminosities com- 
puted from the SED fitting for our sample, we consider Lb 0 i, sed 
throughout the paper but compare the results to those obtained 
with X-ray based Lb 0 i estimates on a larger sample. The Ed- 
dington ratio is then computed as TEdd = Tb 0 i/LEdd, where 
L e dd = 1.26 x 10 38 M B h/A7 o erg s _1 . 


4. Correlation analysis 

In order to test for the dependence of the hard X-ray variability 
on the properties of the selected objects, we correlate the vari- 
ability estimator S v with the X-ray luminosity A 1 4 _ 1 95 ^v, black 
hole mass A7 H h, Eddington ratio /l^id. intrinsic hydrogen column 
density Ah and spectral shape I . We exclude from the correla- 
tion analysis the radio loud sources in order to have a uniform, 
local Seyfert sample for which it is not necessary to correct vari- 
ability and luminosity for redshift and beaming effects, and for 
which the hard X-ray emission is dominated by accretion pro- 
cesses, with no jet contribution. In addition, the results for the 
Seyfert sample will be easily comparable to previous studies also 
focusing on radio quiet objects. 
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Fig. 5. Left: Variability estimator Sv versus Eddington ratio for the Seyfert sample. Yellow triangles represent the average of the variability 
estimator in different Eddington ratio bins. Middle: Photon index versus Eddington ratio for the Seyfert sample. Right: Variability estimator S v 
versus photon index for all the Seyfert galaxies in our sample, separated into the different types. Empty black circles indicate those sources with 
10 23 < N h < 4 x 10 23 cm' 2 . 


We do not find any significant correlation of S v with lumi- 
nosity (Fig. [3]) nor with absorption (Fig. 0). Variability appears 
to be marginally anti-correlated with Eddington ratio when us- 
ing Lboi, sed to compute /t^dd, with a Spearman ranking correla- 
tion coefficient of R con = -0.3 corresponding to a probability for 
chance occurrence of P c orr = 0.03 (Fig. [5} left panel). However, 
the correlation disappears when using a larger sample and the X- 
ray based bolometric corrections to estimate /l^dd- When fitting 
the S v versus Eddington ratio relation with a power law model. 


0 . 09 ± 0.02 


we find Xy « ^Edd ±UU ’ which steepens toXy 11 /l Edd 
when excluding the 3 sources at super-Eddington rates. 

Unlike what is observed in the 2-10 keV band, we do not 
detect any anti-correlation between the 14-195 keV variability 
and the black hole mass (R cor r = 0.06, /\<„- r = 0.64; Fig. [6]). This 
is in agreement with what was found by Shimizu & Mushotzkvi 
(12013 ) on a sm alle r sam ple (squar es in Fig. [6]). On the other 
hand. lCaballero-Garcia et all ( 20121 ) reported a trend of decreas- 
ing variability with increasing black hole mass for the 5 bright- 
est Seyfert galaxies at hard X-rays, with 20-50 keV BAT light 
curves rebinned to 2 days. However, when computing /-\. al for the 
2-days binned, 14-195 keV light curves as presented in our work 
for the same 5 sources, and applying a Spearman rank correlation 
test (more robust than the linear/Pearson correlation one; e.g., 
iPress et~aTll2007h . we do not find any anti-correlation between 
variability and black hole mass (R con = -0.1, P co „ = 0.87). 

The lack of a S v - Mbh anti-correlation is most likely due to 
the fact that the BAT survey samples time scales longer than the 
power density spectrum (PDS) break, 7 b. In fact, the X-ray PDS 
of AGN has a characteristic power law shape with index apos = 
-1 for long time scales and crpos = -2 for times scales shorter 
than a source-dep endent break time 7 b (typically of the order 
of hours to weeks:llJttlevll2005 : Areyakj^t^l JgOOSat iMarkowitzl 


2009. l20HilGonzalez-Martm & Vaughanl2012l) . which is found 
to scale with black hole mass and bolometric luminosity follow- 
ing the relation described by iMcHardv et alJ (2 006 ). This im- 
plies a scaling of the variability amplitude with the same AGN 
properties when time scales of the order or shorter than 7 b are 
considered. Above 15 keV the PDS obtained with BAT data on 
time scales of weeks to years show for 29 out of 30 objects 
a single power law with index crpos ~ -0.8 (SM13), consis- 
tent with the 2-10 keV results and providing an upper limit 
of 26 days for the hard X-ray 7 b of these objects. However, 
if one assumed that 7V, is energy independent and can be esti- 
mated using the lMcHardv et al . ( 2006 ) relation, one would find 



Fig. 6. Variability estimator S v versus black hole mass for the Seyfert 
galaxies in the sample. Empty circles highlight AGN with predicted 
break time larger than 3 years, while blue empty squares point out the 
Seyferts from the SM13 sample. The arrow indicates the Mbh upper 
limit for MR 2251-178. Yellow triangles represent the average of the 
variability estimator in different black hole mass bins. 


that for 12 Seyferts in our sample T^pi-ed is predicted to be 
larger than 3 years. This suggests the possibility that for those 
AGN BAT has measured the variability at frequencies above the 
PDS break, and that therefore a scaling of variability with black 
hole mass should be expectecQ. Yet, no correlation is detected 
(empty circles in Fig. [6]). Recently, based on the 0.2-10 keV 
PDS analysis of 104 AGN and on the detection of a break in 
15 of them. lGonzalez-Martm & Vaughanl(l2012l) recomputed the 
7 b - Mbh - T h<1 | relation finding a weaker dependence on L ho | 
than in the IMcHardv et al. (2006) formulation. Therefore, con- 
sidering the large uncertainties on the mass determination and 
bolometric luminosity estimate, and on the 7 b - Mbh ~ T ho i 
relation, it is not possible to derive from the above arguments 
whether or not this relation is also valid at hard X-rays. 


5 For the majority of the objects in the SM13 sample (20 out of 23) 
r B> pred is of the order or lower than 30 days, consistent with no T B de- 
tected by SM13 in the BAT time frequency range. 
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5. Spectral variability 

5.1. Variability, photon index and Eddington ratio 


A significant anti-correlation is found between variability and 
photon index for Seyfert galaxies, in the sense that AGN with 
harder spectra are more variable (R mn - = -0.5, P c(m - = 2 x 1 0 fi , 
Fig. 0right panel). The photon index used here is the result of a 
simple power law fit to the BAT spectrum. Even when excluding 
the objects with Ah > 3 x 1CL 23 cm~ 2 for which absorption might 
start to have a noticeable influence on the hard X-ray spectrum, 
the anti-correlation is still significant ( P con = 6x 10 4 ). 

A correlation is also identified between photon index and 
Eddington ratio (R corr = 0.56, F corr = 1CL 5 ; Fig. [5] middle 
panel), similar to wh at is observed in the 2-1 0 keV band on 
year time scales (tSobolewska & Panadakisl 120091 ). The best fit 
to the data with a power law model results in T oc d ( Tj l22±0 ' 002 
(r oc dg° 39±0 003 when considering only sources with /lead < 1). 
The correlation is significant also when using the X-ray bolomet- 
ric corrections to compute dedd- This relation has been suggested 
to origin from the geometry of the accretion flow, with stronger 
accreting objects having accretion disks whose inner radius is 
closer to the innermost circular stable orbit and therefore can 
cool more efficiently the hot corona, resulting in softer X-ray 
spectra (tSobolewska & Panadakisll2009) . 

The marginal anti-correlation between variability and Ed- 
dington ratio (Fig. 0 left panel) could therefore be induced by 
the combination of the T - d K dd and S \ - T ones (Sv oc T 2 - 5±0 - 3 
for the total Seyfert sample). 


5.2. Variability in the 14-24 and 35-100 keV bands 

In order to compare the variability at different energies, we use 
the 30-day binned light curves in the 14-24 and 35-100 keV 
bands (see Table 0. These energy bands have been selected to 
provide comparable significance for a typical hard X-ray AGN 
spectrum with power law shape and photon index T = 1.8. 

For a sample of 68 AGN (50 Seyferts, 9 radio galaxies, 9 
blazars), we are able to estimate S v for the two selected bands 
(Fig .|7). In general, the variations are well correlated between the 
two bands (R con = 0.6, P con - = 5 x 10 7 ). However, there is an 
average 10% shift of the Seyfert population towards higher am- 
plitude of the variations at lower energies, aside from a very few 
exceptions. A Kolmogorov-Smirnov test indicates a probability 
Pks > 98% that the variations in the two energy bands are differ- 
ent for the 50 Seyferts in this subsample. This is also true when 
separating the sample into the 13 unabsorbed (Ah < 10 22 cm -2 ) 
and 26 absorbed objects: absorbed and unabsorbed sources show 
the same average variability amplitude, with the lower energy 
band presenting larger variability (( S v)unabs,i 4 - 24 kev = 26 + 2%, 
(S V )abs, 14—24 keV = 28 + 1%, (S v)unabs,35-100keV = 17 + 2%, 
(S v)abs, 35 - 100 keV = 16 ± 1%). On the other hand, on average 
the 18 radio loud AGN do not display a significant dependence 
of variability on energy. 

It is important to stress that the 14-24 and 35-100 keV bands 
have been selected specifically to have similar signal-to-noise ra- 
tio (SNR) for a typical AGN spectrum, therefore excluding the 
possibility that the behavior observed could be due to a bias in- 
troduced by the different uncertainties in the two bands. In fact, 
on average the difference between the SNR in the low- and high- 
energy band is 3% for the total sample (-1.4% for the Seyferts 
and 11% for the radio loud sources), with about half of the sam- 
ple having a larger SNR in the high-energy band. 



14-24 keV S v [%] 


Fig. 7. 35-100 keV versus 14-24 keV variability estimator for the AGN 
in our sample. The dashed line indicates where objects with the same 
variability in the two bands would lie. Even though the variations in the 
two bands are well correlated, for the majority of the radio quiet objects 
the variations in the lower-energy band are larger than those at higher 
energies. Some of the most extreme outliers are labeled. 


When testing the possibility that the additional variability 
in the lower energy band could be ascribed to variations of the 
intrinsic absorption, we find that large and frequent variations 
would be needed. In fact, if we consider a simple power law 
spectrum (T = 1.9) with full-coverage absorption, variations 
of the order of AAh ~ 5 - 10 x 10 23 cm~ 2 would be neces- 
sary to produce flux variations of AEi 4 _ 24 keV ~ 8 - 16% and 
AE 35 -lookeV < 0.7%. This would imply a 10% difference be- 
tween the low- and high-energy flux variations within two ob- 
servations and correspond to a A/*j 4 _. iqji^v ~ 3 — 6% contribu- 
tion to the variability in the full band. However, in order to mea- 
sure an increase in S v by 10% given our sampling baseline and 
rate, unrealistic large ((AAh) at least up to ~ 10 25 cm -2 ) and fre- 
quent (on average over 50% of the light curve points) variations 
in the absorbing column would be required, and would imply 
even larger variations in the energy band <J0 keV. Even though 
variable absorbers are indeed common (iRisaliti et al.l2002h . such 
extreme variations are very rarely observed (see e.g., NGC 1365, 
Risali ti et al.ll2007l) . 


5.3. Hardness ratio versus flux 

In order to compare the spectral variability within the single 
sources, for each object we compute the hardness ratio on month 
time scale, defined as: 


HR = 


E 35-100 keV ~ E[4_24keV 
E 35-100 keV + Fi4_24keV 


( 2 ) 


The dependence of HR on time is then fitted with a constant 
function, a x 2 test is applied, and the corresponding variabil- 
ity index V is computed (as defined in Sect. |3j. Based on 
this, we select a subsample of 18 objects for which V > 1.3, 
as those presenting spectral variability. These are mostly the 
brightest objects in the sample, with a few exceptions (e.g. 
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Fig. 8. Hardness ratio versus 14-195 keV flux for the 18 AGN showing spectral variability (13 Seyfert galaxies on the three top rows and 5 radio 
loud AGN on the bottom row). The black dots are from the 1 -month binned light curves, and red squares indicate the average hardness ratio in 
different flux bins. 


MCG-05-23-016, Circinus galaxy). The relatively low fraction 
of sources showing significant spectral variability is in agree- 
ment with the results presented by iMateos et al l ( 2007 ) who 
found spectral variability to be less common than flux variabil- 
ity on long time scales in the 0.2-12 keV band. In Fig. [8] the 
hardness ratios as a function of flux (normalized to the mean 
flux of each object) are presented. For 8 Seyferts there is a clear 
trend of spectral softening when the source brightens (Mrk 3, 
NGC 3516, NGC 4151, NGC 4507, NGC 4945, IC 4329A, 
NGC 5506, AX J1737.4-2907), while for the other radio quiet 
objects no significant correlation undetected. Our results are in 
agreement with those reported bv ICaballero-Garcia et al.l ( 2012 ) 
for NGC 4945, NGC 2110 and IC 4329, while different trends 
are observed for NGC 4388 and NGC 4151. This is possibly due 
to the different energy bands used to compute HR (20-50 and 
50-100 keV in lCaballero-Garcia et alj[2012l) . 


There is no correlation between flux and spectral variations 
for radio loud objects, with the exception of Mrk 421 and Cen A, 
which exhibit a marginal spectral hardening with increasing flux. 
Cen A has been observed to present a harder-when-brighter 
spect rum also during Suzaku observations dFukazawa et al.l 
1201 ll ). This behavior has been interpreted as due to the increased 
contribution in the brighter state of a hard tail produced in the jet. 


6. Discussion 

6. 1 . Variability of the different AGN classes 

Radio loud objects show larger amplitudes of the variations 
compared to Seyfert galaxies, as expected if their hard X-ray 
emission is dominated by the extreme processes taking place 
in a relativistic jet (e.g.. iLichti et alJl2008l) . while hard X-ray 
emission from Seyfert galaxies is instead dominated by thermal 
Comptonization. This is even more evident considering that the 
variability of higher redshift blazars (4 out of 13 objects have 
z > 0.4) will actually increase when taking into account the ef- 
fect of cosmological time dilation. Among radio loud objects, the 
radio galaxies seem to show rather intermediate variability am- 
plitudes between the blazars and the radio quiet class, in agree- 
ment with the idea that the jet and the coronal emissions are 
probably both contributing. This makes it difficult to disentangle 
the relative importance of these components in the hard X-ray 
spectra of radio galaxie s, as for example shown in C entaurus A 
dBeckmann et alfeoi ill. 3C 111 dde J ong et a l.ll2012l) . as well as 
in other radio galaxies dGrandi et al.ll2006i) . 

Among Seyfert galaxies, there are indications for type 1 .5-2 
objects to be slightl y more variable t han typ e 1, as found already 
at hard X-rays in [Beckmann et al. (2007a). A similar behavior 
has been reported bv ISaxton et al. ( 201 ll) when studying long- 
term X-ray variability properties of a sample of more than 1000 
AGN detected by ROSAT and within the XMM-Newton slew sur- 
vey. Following their study, the Seyfert 2 class has the highest 
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fraction of sources with large variability in the 0.2-2 keV range, 
but it is not clear if the observed variability is due to changes 
in the line-of-sight column density, or it is rather intrinsic to 
the cen tral engine. On the other hand, while in Beckman n et alJ 
(!2007a l there was a tentative trend of more absorbed sources 
being more variable, here we do not observe any correlation be- 
tween variability and absorption column density. Therefore, the 
global properties of variability as a function of AGN class point 
to a general agreement with the basic formulation of AGN unifi- 
cation models, in which the differences are driven in first approx- 
imation only by AGN orientation and by the presence or absence 
of a relativistic jet. 

Within our sample, there are 3 NLS 1 galaxies, a class which 
is known to show peculiar characteristics in the soft X-ray do- 
main, like a steep spectral continuum, and strong and fast vari- 
ability. Most of these distinctive properties of NLS Is have been 
interpreted in the context of relatively small black hole masses 
undergoing rapid evolution, powered by higher accretion rates 
than their broad line equivalent (iPeterson et al.l l20Q0h . The av- 
erage hard X-ray photon index of the NLS1 in our sample 
is indeed slightly larger than for the rest of the Seyfert class 
((Dnlsi = 2.2 + 0.1, (Dseyi = 2.01 + 0.0 2), as deriv ed also 
from INTEGRAL observations on 14 NLS Is (iRicci et al.ll201 lbl : 
iPanessa et al.ll20llh . On the other hand, the average variabil- 
ity amplitude is consistent with what is observed for broad line 
Seyfert 1 ((Sv)nlsi = 15 ± 3%). Even though our NLS1 sam- 
ple is too small to draw any firm conclusion, the discrepancy 
between the soft and hard X-rays results could be due to the dif- 
ferent time scales probed in these studies: while NLS Is do show 
strong variability on short time scales, on long time scales above 
the PDS time break, their timing characteristics do not differ sig- 
nificantly from those of the general type-1 Seyfert population. 
This is in agreement with lSaxton et al.l ( 201 I ). who found sub- 
stantially the same long-term variability in narrow and broad line 
Seyfert 1 within their 0.2-2 keV survey study. 

6.2. Long-term variability, luminosity and black hole mass 


Due to the rather long time-scales probed by BAT, the anti- 


correlations of variability with luminosity 

Barr & Mushotzkv 

119861: Green et al. 1993|: [Papadakis 

2004 

) and with black 

hole mass (iPapadakis & McHardy 

1991 Zhou et all 1201 0t 

Kelly et al. 2011, 2013; Ponti et al. 2012[) often observed below 


10 keV on short time scales are not detected in our hard X-ray 
study (Figs. |3]and [6]). This is in agreement with results obtained 
in the soft X-ray band when year time scales are investigated. For 
example. ISaxton et al.l ( 2011 1 did not find any relation between 
long-term variability and luminosity when studying a 0.2-2 keV 
AGN sample, and argued that this is expected since the sam- 
pled time scales are substantially longer than the typical PDS 
break times of AGN. Moreover, iMarkowitz & Edelsonl (120041) 
analysed a sample of 27 AGN with ASCA and RXTE data in 
the 2-12 keV energy range. Even though they do detect anti- 
correlations between variability and luminosity and variability 
and black hole mass on day time scales, the significance of these 
relations decreases when studying longer time scales (a month 
to years) as well as the slope of the linear fit to the data ap- 
proaches zero, indicating that on long time scales there is no anti- 
correlation observed . iZhan d d20 1 1 1) used RXTE/ ASM data, mon- 
itoring 27 AGN for about 14 years in the 1.5-12 keV band. The 
excess variance computed on the 300-day binned light curves 
does not correlate with the black hole mass. The author argues 
that this points to AGN being in a high/soft state rather than in 
a low/hard state, based on an analogy with Galactic black holes 


(GBH) when solely their timing behavior is considered. In fact, 
PDS of GBHs in low/hard state are observed to present a second 
break at even lower time frequencies, and therefore if a similar 
behavior has to be expected from AGN, a mass dependence of 
the excess variance should be observed also on long time scales, 
which is instead not detected. 

Finally, considering that the BAT PDS on long time scales 
appear to have a slope consistent with apos ~ 1 (SM13), S\r 
computed in this time frequency range is proportional to the PDS 
normalization. Therefore, the lack of correlation between Sy 
and Mbh and Lx indicates that the PDS normalization does not 
seem to scale with black hole mass nor luminosity. Moreover, 
the average variability for Seyfert galaxies is (S v)sey = 19% 
with 68% of the objects having S y values in the range ~ 10- 
35%. The total variability range is smaller by a factor of ~ 6 
with respect to the range measured on shorter time scales by 
iPonti et Si (1201 2h in the 2-10 keV energy band. This suggests 
that the BAT S y . and therefore the normalization of the PDS 
on long time scales, could be very similar in all local Seyfert 
galaxies. 


6.3. Spectral variability 


A significant correlation is detected between variability in the 
14-24 and 35-100 keV energy ranges, with Seyfert galaxies be- 
ing on average more variable in the lower energy band (Fig. [Til. 
The detected correlation suggests that the same process is likely 
to be at the origin of the variability in the two bands. Moreover, 
the presence of such large average variations, even in the 35- 
100 keV range, is incompatible with hard X-ray variability be- 
ing entirely due to variations of the hydrogen column density 
along the line of sight. In fact, as an indication, a variation of 
ANh ~ 5 x 1 0 24 cm 2 would imply only a 3% variation of the 
35-100 keV flux. Yet, a limited fraction of the additional vari- 
ability observed at lower (14-24 keV) energies could be due to 
absorption variations. 

Several studies have detected decreasing variability with 
increasing X-ray energies in the 0.2-12 keV band, either 
through excess variance or PDS amplitude (e.g.. iNandrr 

~>n 2004 “ 


Ponti 


2007 


McHardv et al l 20041; iMarkowitz & Edeiso 
Ponti etalJ 1201 2h . and both on short (< 1 day) and long time 


2001 


scales (> 30 days). The same finding extends to hard X-ray en- 
ergies, with the 14-24 keV band showing on average larger vari- 
ations than the 35-100 keV one in Seyfert galaxies. A similar 
result has been obtained bv iChitnis et al.l ( 2009 ) when compar- 
ing the 1.5-12 keV RXTE/ A SM variability to that measure d in 
the 14-195 keV band by BAT. lShimizu & Mushotzkvl ( l2013l) did 


not detect an energy dependence of the PDS parameters for the 
3 Seyferts for which this investigation was possible, but the cor- 
responding excess variance (computed using the PDS parame- 
ters) decreases through the 14-24, 24-50 and 50-150 keV bands 
they studied. If for observations below 10 keV changes of the 
amoun t of the absorbing material or of its io nisation could play 
a role dRisaliti et alJl2002t IPonti et aTl 120 1 2l) . in the BAT band 
such mechanisms are unlikely to be able to fully explain the ob- 
served properties. 

The variations could instead originate directly at the X-ray 
source. For example, if a power law spectrum pivots at ener- 
gies > 50 keV, the 14-24 keV flux would be expected to vary 
more strongly than the 35-100 keV one dGierlinski & Zdziarskil 
l2005h . If we assume a standa rd Comptonization mode l based on 
a two-phase accretion disk (Haardt & Maraschi 1993), the hard 
X-ray va riability could be driven by changes of the seed pho- 
ton flux (iNandra et alJ 120001 : lArevalo et al] l2005h or by phys- 
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Fig. 9. 14-195 keV versus 2-10 keV normalized excess variance crj ms 
for 10 Seyfert galaxies and 3 radio loud AGN (3C 273, 3C390.3, Cen A) 
in our sample for which 2-10 keV PDS with a measured frequency 
break have been reported in the literature. The dashed line indicates 
where objects with the same variability at soft and hard X-rays would 
lie. 


ical variations in the corona (temperature and optical depth; 
iNicastro e t al. 2000: ILubinski et aD2010h . lPelrucci et alJ feOOO) 
point out that in order to have a pivoting point at high energies 
an increase of the cooling is required rather than a decrease of 
the heating, that would instead determine a low-energy pivoting 
point. 

In addition, the superposition of two components, one con- 
stant (e.g., the reflected continuum) and the other one vari- 
able (for example the Comptonized continuum) could further 
contribute to the observed spectral variations (IShih et al.lf2 002 

tHa - 


Markowitz_&Edelson j2004|ij Pontiet al .1 120061 : IPanadak is et al.l 


2009t Sobolewska & Papadakis 20091) . In the 14-24 keV and 


35-100 keV bands the contribution of reflection to the observed 
flux is about 30-40% for a reflection fraction R = 1 and typical 
AGN spectra with F = 1.6 - 2.1, Nn = 1 - 100 x 10 22 cm 2 , 
and a high-energy cut-off at Ec = 100 - 300 keV. For such spec- 
tra, a larger contribution from reflection is measured in the total 
35-100 keV flux compared to the 14-24 keV one, but only by a 
few percents. Therefore, reflection alone cannot fully explain the 
difference observed in the two hard X-ray bands, even though a 
contribution from it is expected. 

The scenario of a pivoting spectrum superposed to a con- 
stant reflection is in general also consistent with a steepen- 
ing of the hard X-ray spectrum when the source brightens. 
This is indeed observed for the majority of the Seyferts in our 
sample that show some spectral variability (Fig. 0 and fol- 
lows an established trend observed in the X-ray emission of 
numerous^ local radio quiet AGN ({M arkowitz & Edelsonl i2004l; 
ISobolewska & Panadakisl l2009t ICaballero-Ga^ [20121) 

and also in higher-redshift objects dPaolillo et al.ll2004l) . Within 
the Comptonization frame, the softer-when-brighter behavior 
can be understood if an increase of the seed photon power il- 
luminating the corona determines a more efficient cooling of the 
hot electrons, with a resulting drop of the coronal temperature, 
which causes the X-ray spectrum to steepen. 

The correlation between Eddington ratio and photon index 
is_ another known property of AGN samples below 10 keV 
(ISobolewska & Panadakisl 120091) and it is confirmed by our 
stud y at hard X-rays_(Fig. [5j middle panel), in agreement with 
iMiddleton et al.1 ( 2008 ) who suggested this relation to explain 


possible differences in hard X-ray spectral slopes for different 
AGN classes. 

A different behavior arises when comparing the hard X- 
ray to the soft X-ray variability. We collected from the litera- 
ture the PDS parameters for 13 AGN in our sample, i.e., those 
objects with measured soft X-ray PDS and with detected fre- 
quency break JGonzalez-Martm & Vaughanl 120121 SM13, and 
references therein). We compare the normalized excess vari- 
ance cr 2 ms in the 14-195 keV band to the 2-10 keV one, com- 
puted by extrapolating the soft X-ray PDS to low frequencies 
in order to cover the same time-scales as the BAT data set, i.e. 
between v m ; n = 6 x 1 0 9 H z and v max = 4 X 1 0 7 Hz. The 
variability amplitude at hard X-rays is found to be equal or 
larger than that in the soft X-ra ys (F ig. [9]). This is consistent 
with what has been found b vICaballero-Garcia et al.l ( 20121 ) and 
IShimizu & Mushotzkvi (1201 31) . yet it is surprising, as it seems to 
be an opposite trend to those observed within the soft and hard 
X-rays bands alone. One reason could be that, when extrapolat- 
ing the 2-10 keV PDS to year time-scales, the uncertainties on 
the measurement of the PDS slope introduce a large scatter in 
the estimated excess variance. Indeed, for 6 out of 13 objects 
the low-frequency PDS slope has been fixed to a = -1 (from 
iGonzalez-Martfn & Vaughanll2012|) . Slightly steeper slopes by 
Arr = 0.1- 0.2 would imply a comparable or larger variability at 
soft rather than at hard X-rays. 

If instead the hard X-rays are indeed more variable than the 
soft X-rays on long time scales, this might indicate an even more 
complex dependence of variability with energy. The difference 
between soft and hard X-ray variability could be, for example, 
related to the size of the emitting region. If the high-energy 
plas ma is located inside the accretion disk or above its inner 
part (IZdziarski et alj[l999t ILubinski et al.ll20lol) . while the low- 
energy emitting region is associated to some outflow from the 
disk or a second plasma with lower temperature dPetrucci et al.l 
1201 3l) . then the high-energy emission could be more variable 
than the lower-energy emission. However, the change in relative 
dominance of the two components would have to happen exactly 
at energies of 10-15 keV. In addition, it would be necessary to 
explain why within the soft and hard X-ray bands this trend is 
inverted, with the higher-energy emission being in general less 
variable than the lower-energy emission. 


6.3.1 . Variability and photon index anti-correlation 

Another interesting difference compared to the lower energies 
and lower time frequencies is the trend of harder sources being 
more variable at hard X-rays (Fig. [5] right panel). In fact, several 
papers report the opposite trend of softer sources being more 
variable on time scales of the order or below a day, even though 
with different degrees of confidence. lO’Neill et al.l (120051) found 
a marginal correlation betwee n ph oton index and exc ess vari- 
ance, while iGreen et all (1 19931) and iGrune et~ai] (1200 ll) found a 
more significant correlation in the 1-8 keV ( EXOSAT ) and 0.2- 
2 keV band ( ROSAT ), respectively. One has to keep in mind 
though that a contribution from absorption variations could be 
important in the ROSAT data. The largest soft X-ray samp le 
for this kind of study has been presented by IPonti et al . ( 201 2). 
They confirm the correlation between excess variance and pho- 
ton index, which they ascribe to the the fact that accretion rate 
correlates with both variability and photon index in their sam- 
ple. On the contrary, at hard X-rays and low time frequency, it 
seems more likely that the correlation between accretion rate and 
photon index and the anti-correlation between photon index and 
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Fig. 10. Measured photon index versus intrinsic absorption for the 
Seyfert galaxy sample. The lines correspond to the measured photon 
index obtained from a simple power law fit of simulated power law 
spectra with intrinsic photon index T = 1.6, 1.9 and 2.2 (from bottom 
to top) modified by the corresponding Ah- Upper limits on Ah are indi- 
cated with arrows. 


variability are at the origin of the marginal anti-correlation be- 
tween accretion rate and variability (Fig. 0). 

The photon index versus variability anti-correlation is in- 
stead commonly observed in Galactic black holes (GBH), which 
show larger variability during the hard state, i.e., when the 
corona rather than the accretion d isk emission is believed t o 
dominate the X-ray spectrum (e.g.. iMunoz-Darias et al.ll201 ll) . 
Moreover, a decreasing variability with the softening of the spec- 
trum is detected in the so-called hard-intermediate state, when 
the less variable accretion disk component starts contributing to 
the X-ray spectrum of GBHs dBelloni et al.l201 ll) . This might in- 
dicate that the larger variability in the harder BAT sources could 
be associated to their stronger Comptonization. 

The anti-correlation between variability and photon index 
in our sample could be, at least partly, induced by the effect 
of different reflection fractions. On one hand, if one assumes 
that the bulk of reflection is mostly constant, objects with larger 
reflection would be expected to be less variable. Constant or 
slowly variable reflection components observed in a number of 
objects have been explained either due to a distant reflector, 
with l ight travel time dam ping out variability (e.g. JMiller et alJ 
i 2008l ' Bianchi et alii 2009 ). or as a consequence of strong gravi- 
tational light bending in the close vicinity of the black hole (e.g ., 
iMiniutti & Fabianll2004 iPonti et aJ1l2006l iMiniutti et al.l l2007). 
On the other hand, there are indications that objects with larger 
reflect ion have softer X-ray spectra. In fact, Zdziarski et al. 
( 1999 ) reported a strong correlation indicating that intrinsically 
softer sources have larger reflection fractions, possibly due to 
the reflector being also the source of seed photons for the 
Comptonization emission. In addition to this effect, when fitting 
with a simple power law model spectra with different reflection 
strengths and the same intrinsic continuum, the measured pho- 
ton index in the 14-195 keV band results steeper for larger re- 
flection, with AF <7% (depending on the intrinsic continuum 
shape and on the position of the high-energy spectral turnover) 
for a reflection fraction R changing from 0 to 2. Reflection frac- 
tions R > 1 might be the result of peculiar geometries (causing, 
e.g., high covering factors; Malzac 20 0T1). of light bending ef - 
fects in the vicinity of the black hole^Miniutti & Fabian 2004), 


or of a variable continuum combined with a delayed answer 
of the reflected compone nt, due to the distance of the reflector 
dMalzac & Petruccill2002l) . 

When the accretion disk is considered as the main reflec- 
tor, a larger reflection component is expected from face-on com- 
pared to edge-on AGN. Indeed, when dividing our sample into 
the different AGN classes, Seyfert 1 appear at the higher end 
of the photon index distribution, consistent with larger reflec- 
tion b eing asso ciated to softer spectr a. Furthermore. iRicci et alJ 
(1201 lhl) and IVasudevan et alJ (1201 3h have recently shown that 
objects with absorption in the range 10 23 < Ah < 10 24 cm~ 2 
present a stronger reflection component than those at lower Ah- 
When looking at the distribution of photon indexes as a func- 
tion of Ah, sources in the range 10 23 < Ah < 4 x 10 23 cm 2 
have a softer spectrum than the rest of the Seyfert 2 population 
and comparable to the average Seyfert 1 one (Fig. ITOl) . possi- 
bly pointing out a stronger reflection component. We note that 
the photon index drop at Ah > 4 x 10 23 cm~ 2 cannot be simply 
explained by the fact that our fitting does not take into account 
absorption. In fact, when simulating absorbed spectra and fitting 
them with a simple power law, the measured photon index de- 
creases much less dramatically as a function of Ah then observed 
in the data (see lines in Fig.lTOl). 

However, caution should be used considering that for highly 
absorbed sources simple absorption models fail to correctly re- 
produce the spectra in presen ce of large reprocessed emissio n 
(Lubinski et al. in preparation. lYaqoobl2012llBraito et alj2013i) . 
and that stacked spectra might be influenced by the effects of 
combining objects with different spectral shapes contributing 
with varying weight to the different energy bands. In addition, 
the validity of the T — R correlation reported by IZdziarski et al.l 
is still debated. Some later studies confirmed this re- 


lation (Mattson et al. 2007|_ Dadina 

2008; Winter et al. 

2012b 

whereas others did not (iMolina et al. 

2009; de Rosa et al. 

20 1J), 


pointing out the difficulty in disentangling a real physical rela- 
tion from the degeneracies between photon index, reflection and 
high-energy cut-off in the commonly us ed X-ray fitting models 
dMattson et al.ll2007t I Winter et~akll20 1 2l) . 

An alternative explanat ion could come from a sim ilar trend 
observed in NGC 7469 (iNandra & Panadakisl 1200 ll) . which 
presented a highe r excess variance when the spectrum of the 
source was harder. iNandra & Panadakis ; ( 200 1) suggest that this 
and other timing properties of this AGN could be explained if 
the inner regions of the corona are hotter than the outer ones. 
The hard X-ray emission would then be produced in the inner- 
most parts of the corona and therefore, be more variable and 
on shorter time scales than the soft X-ray one, produced in the 
larger and cooler outer zone. However, these results are based 
on a 30 days long observation at <15 keV, hence the time scales 
and energies are different than those of our work. Similarly, 
in the model first pro posed by Lyubarskiil ( 1997b and then 
recon sidered by others (IChurazov et al.ll200lb lArevalo & Uttlevl 
2006), different variability patterns at different energies and 
on different time scales are explained with accretion disk 
fluctuations propagating inward and illuminating a corona with 
emissivity gradient, which would therefore modulate the X-ray 
source. 


7. Conclusions 

We have characterized the 14-195 keV variability of the largest 
sample of hard X-ray selected AGN provided by Swift/BAT up 
to date. The amplitude of the variations has been quantified and 
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their dependence on energy and on basic AGN properties, such 
as black hole mass, luminosity, Eddington ratio and intrinsic ab- 
sorption, has been investigated. The average characteristics of 
hard X-ray variability are in general agreement with the uni- 
fied models, indicating larger variability in jetted sources and no 
clear dependence of the variations on AGN orientation. The am- 
plitude of the hard X-ray variations and its spectral properties in- 
dicate that variability cannot be driven by absorption variations. 
On the other hand for radio quiet AGN it is rather consistent with 
changes of a power law continuum over a constant component, 
with the power law pivoting around > 50keV. Such variations 
of the continuum emission would be achieved for example with 
a change in the corona temperature, determining a shift of the 
high-energy cut-off at hard X-rays, either due to intrinsic coro- 
nal variations or to a change in the seed photon flux. 

In general, variability at hard X-rays has very similar proper- 
ties to that in the softer X-ray band, when similar time scales are 
compared, compatible with variability being driven by the same 
mechanism across the full X-ray band. However, some discrep- 
ancies, like AGN with harder spectra being more variable and 
the apparent lower variability in the 2-10 keV band, could be as- 
cribed to additional components in either band (e.g., reflection, 
absorption), and in any case could provide interesting insights 
on the intrinsic differences of the emission in the soft and hard 
X-ray bands, as for example concerning the emitting regions and 
their sizes. 

The lack of dependence of the long-term variability on black 
hole mass and luminosity, and the relatively narrow range of 
variability measured in Seyfert galaxies suggest that the normal- 
ization of the X-ray power spectrum is independent on black hole 
mass and luminosity and might have a similar value for all local 
Seyfert galaxies. 

Thanks to the SwiftfBAT on-going observations it will be 
possible to further extend the studied time scales and increase 
the available AGN sample, obtaining the best hard X-ray AGN 
variability sample for some time to go. At the same time, 
pointed observations with Suzaku, INTEGRAL and NuSTAR will 
continue to investigate the spectral variability in bright and 
well studie d AGN as already successfully done in a number 
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Table 1. Average properties of the BAT AGN in our variability sample, as a function of the AGN class and the amount of intrinsic absorption. 


Type 

(z) a 

(F 14-195 kev)“ 

<ry 

(Sy) a 

(S v,14-24kev) fc 

(S V, 35-100 kev)* 



[1CT 11 ergcm~ 2 s _1 ] 


[%] 

[%] 

[%] 


Sey 1 

0.033 (36) 

6.4 ± 0.8 

2.01 ±0.02 

16.5 + 0.9 

21 ±2(15) 

13 + 2 

Sey 1.5 

0.017(17) 

10 + 2 

1.92 ±0.03 

24+ 1 

32 + 2(11) 

23 + 2 

Sey 2 

0.016(32) 

10+ 1 

1.86 ±0.03 

19.9 + 0.7 

28 ± 1 (22) 

16.2 + 0.9 

NLS1 

0.019 (3) 

11+6 

2.2 ±0.1 

15 + 3 

24 + 4 (2) 

9 + 4 


unabs Sey 
abs Sey 
CT Sey 

0.031 (41) 
0.016 (46) 
0.011 (6) 

6.7 ± 0.7 
10+ 1 
13+4 

2.00 ±0.02 
1.90 ±0.03 
1.8 ±0.1 

18.5 + 0.8 
19.7 + 0.6 
16 ± 1 

25 + 2 (23) 
28 ± 1 (26) 
37 + 6 (4) 

17 + 2 
16.4 + 0.8 
9 ± 1 

all Sey 

radio galaxies 
blazars/quasars 

0.023 (88) 
0.033 (9) 
0.855 (13) 

8.5 ±0.8 
22 ± 12 
7 + 3 

1.95 ±0.02 
2.1 ±0.2 
2.0 ±0.2 

19.3 + 0.5 
24+ 1 
33 + 2 

26.9 ± 0.9 (50) 
26 + 2 (9) 

36 + 2 (9) 

16.7 + 0.7 
26 + 2 
35 + 2 


Notes. In parentheses, the number of objects used for the given average value is indicated. “ All objects in the sample have measured redshift and 
photon index (from a simple power law fit). b Only objects having variability measured in both the 14-24 and the 35-100 keV bands have been 
considered here. 
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Table 2. Variability estimator S v > black hole mass and bolometric luminosity for the selected SwiftfBAT AGN sample. 


Source name 

Type 

X 

S y 

S V, 14-24 keV 

5 V, 35-100keV 

log M B h 

log Z^bol 



[ 10~ 4 c/s/det] 

[%] 

[%] 

[%] 

[l^o] 

[erg s' 1 ] 

QSO B0033+595 

BLLac 

0.60 

99+19 

94+ 16 

127 !j 2 

io!f 



NGC 235A 

Seyl 

0.74 

14^ 


8.8 ±0.5° 


Mrk 348 

Sey2 

2.56 

25 + 3 

30 + 4 

22+4 

7.2 ± 0.7® 

46.4-1 

Fairall 9 

Seyl 

0.84 

io! 4 



8.9 ±0.5" 

44.9° 

NGC 526A 

Seyl. 5 

0.91 

35+6 

43 - + 10 

38:;® 

26 + 10 

8.0 ± 0.5" 


NGC 612 

NLRG 

0.76 

27^ 

9 + 5 

8.5 ±0.5" 


ESO 297-018 

Sey2 

1.01 

16+4 

27 + 8 


9.7 ±0.5* 

44.21 

NGC 788 

Sey2 

1.19 

15+4 

19 -a 


7.5 ±0.7® 

44.3-1 

LEDA 138501 

Seyl 

0.88 

24 + 7 

17+8 

1 ' -9 

24+12 

Z ++13 



NGC 931 

Seyl. 5 

1.07 

12 + 5 


7.6 ±0.3® 

44.5-1 

[HB89] 0241+622 

Seyl 

1.38 

15 -5 

18 + 8 

14!® 



NGC 1 142 

Sey2 

1.31 

30 + 4 

31 

30 !® 

9.4 ±0.5® 

44.81 

NGC 1275 

RG 

1.65 

15*5 

14 + 4 

25 + 13 

^ -14 

8.5 ±0.7® 

45. 0l' 

NGC 1365 

Seyl . 8 

1.09 

22 + 4 

26*® 

10 + 5 

7.7 ±0.3® 

43.81 

ESO 548-G081 

Seyl 

0.68 

36 + g 

51 +16 
J -15 

33 ± 14 

8.9 ±0.5" 

44.01 

3C 111.0 

RG 

1.91 

16 + 4 

15 + 6 

20 *o 

21 * ? 
Z,A -8 

9.6 ± 0.8® 


3C 120 

RG 

1.56 

23+4 

23*® 

7.83 ± 0.04 c 

45.3-1 

UGC 03142 

2MASX J04440903 +28 13003 

Seyl 

Sey2 

0.77 

0.81 

34 + 9 
J -10 

91 ± 17 

56 !£ 
104*^0 
20 *;, 

73+2i 

' -20 



1RXS J045205. 0+493248 

Seyl 

1.03 

21 

8.0 ±0.3"' 

45.0® 

2MASX J05054575 -2351139 

Sey2 

0.87 

33+6 

28 * 8 

16 + 7 

7.5 ± 0.5" 


IRAS 05078+1626 

Seyl. 5 

1.50 

17 1 

3 + 2 


/r qq + 0.01 c 

6.98 _ 02 


Ark 120 

Seyl 

1.07 

9 -5 



0 1 q +0.05 ^ 
8,18 -0.06 

44.9-1 

ESO 362-18 

Seyl. 5 

0.76 

9+4 

19 + 9 


9.0 ± 0.5" 


PICTOR A 

BLRG 

0.63 

5 1 *9° 

29 

^ -10 

33 +14 

J -13 

69 +18 

u -17 

7.6 ± 0.5" 


PKS 0521-36 

BLLac 

0.52 

55 ± 17 


8.7 ± 0.37 


PKS 0548-322 

BLLac 

0.55 

24 +9 
^ -10 

15 :i 

30 +15 

-14 

8.2 ± 0.3-7 

44.1-1 

NGC 2110 

Sey2 

4.77 

33 + 3 

35+4 

32 + 4 

8.3 ±0.3® 

MCG +08-11-011 

Seyl. 5 

2.15 

33 + A 
29 

-10 

32 + 5 

31+6 

8.1 ± 0 . 6 ® 

45.01 

2MASX J05580206-3820043 

Seyl 

0.60 

35* u 

-12 


8.4 ±0.5" 


IRAS 05589+2828 

Seyl 

1.14 

33 + 8 

25 ± 10 



44. 7 d 

Mrk 3 

Sey2 

1.86 

35 + y 

54+ 20 

J -15 

29 !® 

8.7 ±0.3® 

44.5-1 

Mrk 6 

Seyl. 5 

0.93 

18 1 


12 !® 

8.13+0.04" 

44.31 

Mrk 79 

Seyl. 2 

0.75 

6+4 

29 +9 
-10 


8.4 ± 0.5" 

44.6-1 

2MASS J075941 8 1 -3843560 

Seyl. 2 

0.82 

10 + 5 

31 ± 12 

17+8 

8.3 ±0.5® 


Mrk 1210 

Sey2 

0.91 

24*g 

9 + 6 

40+14 

6.5 ± 0.7 s 


Fairall 272 

Sey2 

0.64 

50 + * 2 

JU-13 

47!“ 

36 !n 



[HB89] 0836+710 

blazar 

0.90 

12 + 5 

12 !® 

22 i 

9.4 ±0.3® 


MCG -01-24-012 

Sey2 

0.70 

19+9 
1 -10 

26 -n 


7.2 ± 0.5" 


MCG +04-22-042 

Seyl. 2 

0.63 

27 * 8 

-li 

37 + 12 

-14 

8.5 ±0.5" 

44. 5 d 

Mrk 110 

NLS1 

0.94 

12 + 5 

24 + 9 


7 42+0.09* 

- 0.1 , 
6.3 ±0.5® 

44.7-1 

MCG -05-23-016 

Sey2 

3.44 

19 + 3 

23+4 

16 + 5 


NGC 3081 

Sey2 

1.28 

19*® 

I 8-10 

17+9 

7.4 ±0.3® 


NGC 3227 

Seyl. 5 

1.81 

27 * 5 

> 
y -9 

33 -a 

25+6 

6.9 ±0.1' 

43.9-1 

NGC 3281 

Sey2 

1.37 

21 i 

3+2 

8.0 ±0.5® 

44.21 

Mrk 421 

BLLac 

3.17 

88 !§ 

107 ± 12 

8.3 ±0.3® 


NGC 3516 

Seyl. 5 

1.77 

26+4 

39 + 5 

20! 4 

7 50 +0 04? 
1 -0.06 

44.3-1 

NGC 3783 

Seyl 

2.84 

15*3 

21 i 

20 + 3 

9 + 4 

7 47 +0 07* 
' 1 -0.09 

44 . 4 I 

UGC 06728 

Seyl. 2 

0.51 

15 1 


6.4 ±0.3"' 

44.3 rf 

2MASXJ1 1454045-1827149 

Seyl 

0.86 

24! 
24 *g 

25 ± 13 


6.7 ±0.5" 


LEDA 38038 

Sey2 

0.85 

30 +12 
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Table 2. continued. 


Source name 

Type 


Sy 

S V, 14-24 keV 

S V, 35-100keV 

log M B h 

log L\) oi 



[ 10~ 4 c/s/det] 

[%] 

[%] 

[%] 

[Me] 

[erg s' 1 ] 

NGC4151 

Seyl.5 

8.06 

29 3 

30 + 4 

28 + 3 

7 5+°.!* 

43.7* 

NGC 4388 

Sey2 

4.01 

26 + 2 

26! 3 

25 + 3 

7.2 ±0.6* 

43.7* 

3C 273 

quasar 

6.23 

29 + 3 

29 + 3 

28 + 3 

9 81 +(U b 

y '° l -0.07 

47.1° 

NGC 4507 

Sey2 

2.85 

20 + 4 

30 + 5 

15 +3 

4 

7.6 ± 0.6* 

44.3* 

ESO 506-G027 

Sey2 

1.28 

43+7 

43 ± 10 

31+7 

8.6 ±0.5* 

44.6* 

NGC 4593 

Seyl 

1.40 

15+4 

17!® 


6.99 + _Z 8b 

44.1*' 

WKK 1263 

Seyl.5 

0.72 

36 + _l 
36%° 

54+19 

40 +14 

-15 

19 
y 10 

8.0 ± 0.5* 

44.3* 

SBS 1301+540 

Seyl 

0.52 


7.5 ±0.5“ 

44.9* 

NGC 4945 

Sey2 

3.48 

34 + 4 

48 ± 10 

28+4 

6.2 ±0.3* 


ESO 323-077 

Seyl. 2 

0.73 

19 + g 

42 +20 
-21 


7.4 ± 0.6* 


Cen A 

RG 

20.40 

39+4 

37 + 4 

39 + 4 

8.0 ± 0 . 6 * 


MCG -06-30-015 

Seyl. 2 

1.26 

13 

21 ^9 


r t+ 0 . 1 & 
6 - 7 - 0.2 

44.1* 

NGC 5252 

Sey2 

1.65 

59!® 

61%° 

54 + 8 

9 03 +0 ' 4 
- 0.02 

45.4*' 

4U 1344-60 

Seyl.5 

1.67 

14 + 4 

23 + 7 

7 -5 

9_ + 

7.4 ±0.1* 

44.1* 

IC 4329A 

Seyl. 2 

4.82 

17 + 3 

31+5 

8.1 ± 0 . 2 * 

44.8* 

Mrk 279 

Seyl.5 

0.63 

36!$ 

o 'j + 12 

^ ' -13 

35 +u 
J -12 

8.6 ±0.5“ 

45.0° 

Circinus Galaxy 

Sey2 

4.54 

11 ± 1 

17 + 2 

6.3 ±0.1' 


NGC 5506 

NLS1 

4.12 

15 +3 

31 * 
JL -4 

25 ! 4 

19 +' 
-8 

8+4 

6.7 ± 0.7* 

44.1* 

NGC 5548 

Seyl.5 

1.19 

24 !® 

7.82 ± 0.02* 

44.8* 

ESO 511-G030 

Seyl 

0.76 

33 + 11 

53 +15 
J -16 

39 +lS 

-24 

25 +l2 

89^ 

„+]$ 

JZ -ii 


8.7 ± 0.5* 

44.4* 

NGC 5728 

Sey2 

1.35 

27 !® 

19!® 

8.5 ± 0.5“ 

43.7* 

Mrk 841 

Seyl 

0.64 

21 ± 10 


8.5 ±0.7* 

45.8*' 

PKS 1510-08 

quasar 

0.83 

41 !j? 

27 

8.7 ± 0.5* 

46.4* 

VII Zw 653 

Seyl. 2 

0.87 

21 

24+13 

^-16 



Mrk 1498 

Seyl. 9 

0.74 

5 + 3 


8.6 ± 0.5“ 

45.3' 

2MASX J 1648 1523-3035037 
NGC 6240 

Seyl 

Sey2 

0.77 

1.06 

47 + 8 
22 ^® 

54 +16 

-2Q 

54 +1 ' 
J -18 

42 +15 
^-16 

4 + 3 



Mrk 501 

BLLac 

0.89 

61+8 

64 + 9 

sr\ +12 
62 _ 

21 +iO 

-A? 

35!“ 

9.2 ±0.3* 


1RXS J165605. 6-520345 

Seyl. 2 

0.69 

21 ± 10 

r\sr +15 

44^ 

^-22 

7.9 ±0.5* 


2MASS J 1656 1677-3302 127 

blazar 

0.87 

41+7 



NGC 6300 

Sey2 

1.65 

17 + 3 

18^ 

8! 4 

5.5 ± 0.4* 

42.9* 

AX J1737.4-2907 

Seyl 

1.96 

15 +3 

1 J -4 

33 + 7 


8.9 ±0.7* 


IC 4709 

Sey2 

0.66 

r\r\ +10 
ZZ -11 

34-Ig 




PKS 1830-21 

FSRQ 

1.04 

18 + g 

■it! 

38-21 




3C382 

Seyl 

1.47 

15 +5 

11+5 

9.2 ±0.5* 

45.6° 

ESO 103-035 

Sey2 

1.88 

12 + 3 

17 + 5 


7.1 ±0.6* 

44.6* 

3C 390.3 

BLRG 

1.69 

17 + 3 

16 + 5 

12 + 5 

o AC.+ 0.09 b 

8.46 _ 0 j 

44.9* 

Fairall 5 1 

Seyl 

0.70 

20 + 8 
2*3 

45 +17 

H -19 


7.5 ± 0.5 m 


NGC 6814 

Seyl.5 

1.19 

oz:+14 

36_i7 

i9!; 

7.1 ±0.2* 

43.9* 

Cygnus A 

RG 

2.13 

31 +4 

35^ 
J -in 

20 ^ 

32 +5 

35+6 

9.4 ±0.1* 

45.7* 

QSO B1959+650 

BLLac 

0.69 

42 + 9 


8.1 ±0.3* 


MCG +04-48-002 

Sey2 

1.13 

24 + _l 




4C +74.26 

Seyl 

0.87 

21 !® 

141 

28 + 11 

9.6 ± 0.5* 

46.2* 

Mrk 509 

Seyl. 2 

1.60 

17 + 7 

16!® 

12! 7 6 

o 1 r +0.03 * 
8 - 16 -0.04 

45.0* 

IC 5063 

Sey2 

1.07 

18 + 5 

38 ± 13 


7.7 ±0.5" 

44.5* 

2MASX J21 140128+8204483 

Seyl 

0.55 

12 + 6 

o r +15 
36 -L 6 

29 ! 3 


8.8 ±0.5* 

45.5* 

28P 206 

RG 

3.20 

30 + 4 

31+5 



SWIFT J2 12745. 6+565636 

NLS1 

0.73 

20^8 

21+9 

20+11 

7.2 ± 0.5“ 


6dF J2 132022-334254 

Seyl 

0.72 

40 + 7 

49 + 12 

H -13 

0/7 +15 
36 16 



PKS 2149-306 

FSRQ 

1.07 

26 + 7 

29 + 10 



NGC 7172 

Sey2 

2.63 

28 + 3 

33 + 7 

29 + 5 

7.7 ±0.6* 

43.8* 
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Table 2. continued. 


Source name 

Type 

X 

S\ 

S V, 14-24 keV 

S V, 35-100keV 

log M B h 

log Lbo\ 



[10~ 4 c/s/det] 

m 

[%] 

[%] 

[M q ] 

[erg s' 1 ] 

NGC 7213 

Seyl.5 

0.68 

23 +9 
^ -10 



8.6 ±0.5" 

44.3-^ 

NGC 7314 

Seyl 

0.80 

21 +9 
-10 



6.0 ± 0.5* 

43. 0 7 

3C 454.3 

FSRQ 

1.74 

59 + 8 

6o 
11 s 

56 + 9 

9.2 ±0.7* 

47.3 / 

MR 2251-178 

Seyl 

1.67 

lO! 4 

1 ± 1 

<6.9* 

45.8° 

NGC 7469 

Seyl 

1.10 

23+6 

22 i, 

20!|, 

7.09 ± 0.05* 

45.3 / 

Mrk 926 

Seyl.5 

1.81 

27 : 4 

26 + 7 

21 +5 
Z,A -6 

7.1 ±0.6* 


NGC 7582 

Sey2 

1.21 

23+5 

24 :i 

12 + 6 

8.3 ±0.5" 

43. 3 ; 

NGC 7603 

Seyl.5 

0.82 

U+6 

-7 



8.1 ±0.3* 

44.7 / 


Notes. The variability estimator is computed in the 14—195 keV (Sv), 14-24 keV (Sv, i 4 - 24 kev). and 35-100 keV (5 v, 35 -iookev) bands, x is the 
average count ra te in the full 14-195 k eV band. The bolometri c luminosities have been estimated from the fi tting of the spectral ener gy distribution 
see S ect . [331. "IWinter et all d2009h: ^ [B eckmann et all d2 009h and refe re nces therein; c [Grier et al]d2012|); "IVasudevan et al.ld2009|); e IStalin et al.1 
20lTh: /|WoS & Un\dd200f); 1 'lBian & Gid d2007l); '’iKaspi et alJd2007h. ' iDennevet all <20 1 Cl): 'IVasudevan et alJd2010h : * Nikolai uket alJ <20041) : 
iGreenhill et al.l d2003h r™ lPadovani & Rafanellil ll988): " Mali ziaetakl <2008h : " IVasudevan & Fabianl d2007li . 
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